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Abstract 

Superradiant instability turns rotating astrophysical black holes into unique probes of 
light axions. We consider what happens when a light axion is coupled to a strongly coupled 
hidden gauge sector. In this case superradiance results in an adiabatic increase of a hidden 
sector CP-violating ^-parameter in a near horizon region. This may trigger a first order 
<^ phase transition in the gauge sector. As a result a significant fraction of a black hole mass 

is released as a cloud of hidden mesons and can be later converted into electromagnetic 
^ radiation. This results in a violent electromagnetic burst. The characteristic frequency of 

Q such bursts may range from ~ 100 eV to ~ 100 MeV. 

m 

1 Introduction 

p I String theory landscape points towards ultimate unification of particle physics and geography. 

If the landscape is real it is very likely that in the short and midterm future particle physicists 
^ will be busy by mapping the compactification manifold where we happened to live. In the most 
£j optimistic scenario — if extra dimensions are large — this can be achieved directly at the LHC and 
^ future colliders. However, even if extra dimensions are small they are still likely to be populated 
j> by rich structures accessible for collider, laboratory, astrophysical and cosmological probes at 
^ relatively low energies. These structures escaped detection so far as a consequence of their distant 
QQ geographic location in extra dimensions. 

The structures which are natural to expect in the landscape on theoretical grounds include 
CN hidden gauge sectors ("hidden valleys", [1]) and a plenitude of axion-like particles, dark photons 
^ and photini (the "axiverse" , [2] ) • Confinement scales of strongly coupled gauge sectors and axion 
^ masses are exponentially sensitive to the parameters of the compactification manifold and may 
^ take arbitrarily low values. For axions there is also an observational indication that they may 
^ be light. Indeed, not only the QCD contribution to the mass of the QCD axion is very small 
5-j (~ 10~^° eV for the GUT scale axion decay constant), but also additional string corrections to 
the axion potential, expected to be there on general grounds, have to be suppressed at least by 
extra ten orders of magnitude. Consequently, one may expect the presence of even lighter axions 
with masses dominated by these corrections. 

A lot of efforts were concentrated recently on exploring the possibilities to discover these struc- 
tures in the near future experiments and observations (see, e.g., P-[9]). The focus of the current 
work will be the intriguing possibility to probe light axion particles with masses in the range 
10"^'' -r- 10"^'' eV with ongoing observations of astrophysical black holes [21 [10]. This possibility 
is related to the famous Penrose process allowing extraction of kinetic energy from rotating black 
holes If a boson with a Compton wavelength of order the size of a black hole is present. 
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the Penrose process results in a superradiant instability [I6]-[T9]. A black hole releases its ro- 
tational energy by creating a cloud of a rotating Bose-Einstein condensate in the near horizon 
region. Compton wavelengths of axions in the above mass range match sizes of astrophysical black 
holes and Ref. [10] studied opportunities provided by the superradiant instability for the discovery 
of these particless. The focus of Ref. [10] was mainly on model-independent purely gravitational 
signatures of superradiance. Perhaps one of the most interesting conclusions is that Advanced 
LIGO may be able to discover the QCD axion if the decay constant is high. 

Instead, here we will consider model-dependent signatures arising when an axion causing the 
instability is coupled to a QCD-hke hidden gauge sector (a hidden valley). The presence of these 
signatures depends on some of the details of the matter spectrum in the hidden valley. The payoff 
that we get for introducing a certain amount of model-dependence is that the resulting signals are 
electromagnetic rather than purely gravitational. 

Concretely, we are exploiting the following effect. From the point of view of a hidden gauge 
sector a growing cloud of axions surrounding a rotating black hole acts as a QCD laboratory with a 
slowly oscillating CP-violating ^-parameter. The amplitude of oscillations grows and reaches order 
one at the late stages of the instability. As we discuss in section [2} it is a very common feature of 
gauge theory dynamics that under an adiabatic increase of 9 the vacuum which is stable at ^ = 
first turns into an overheated metastable minimum and later on either becomes perturbatively 
unstable or even totally disappears as an extremum of the potential. 

Unfortunately, this possibility is not realized in the real-world QCD. However, it may be 
exhibited in one of hidden valleys. If this happens, the superradiant instability will trigger an 
avalanche phase transition in the hidden valley and a cloud of strongly interacting mesons will 
get produced. The total energy of the cloud depends on the axion decay constant, and for high 
scale axions may constitute 10~^ -r- 10~^ of a black hole mass. Unlike axions, whose couplings to 
the rest of the matter are strongly suppressed by the decay constant, hidden sector mesons may 
have substantial interactions with the Standard Model fields. The net result is that a significant 
fraction of a black hole mass can be converted into electromagnetic energy on a characteristic time- 
scale set by a black hole size. This would give rise to violent transient astronomical sources — 
electromagnetic echoes of hidden valley avalanches — that can be looked for in a wide range of 
frequencies. 

The rest of the paper is organized as follows. In section |2] we review how a vacuum structure 
with a non-trivial dependence on a 6'-parameter may arise in a QCD-like theory. As a fully 
calculable example we use a QCD-like theory with light quark flavors, where the relevant 
dynamics can be analyzed with the help of a low energy chiral Lagrangian. In section |3] we discuss 
in details how superradiant instability may trigger an avalanche phase transition in a hidden valley. 
In section |4] we estimate characteristic energies, frequencies and time-scales of the corresponding 
electromagnetic signal as a function of an axion mass. In the concluding section [5] we present our 
conclusions and outline directions for the future work. 
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2 Phase Structure of Hidden Valleys 



To illustrate how an adiabatic change of a 6'-parameter can trigger a phase transition in a strongly 
coupled gauge sector, let us consider SU (Nc) QCD-like gauge theory with light quark flavors 
(this section is mostly a summary of well-known facts, an incomplete list of original references 
include |20l |2T1 122] ) . In this case all the relevant dynamics can be analyzed in a weak coupling 
regime using the low energy chiral Lagrangian. Indeed, from the chiral perturbation theory we 
get the following effective Lagrangian 

A3 

C = ^Tr dU^dU + —Tr (Me-^^/^f/ + Me'^^^U^) , (1) 

where U G SU (N) are pions, M is the quarks mass matrix and the pion decay constant is of 
the same order as the QCD scale A. With appropriate field redefinitions one can always make M 
real and diagonal, M = diag{mi,m2 . . . "hitv}, with 

mi > m2 > . . . wijv > . (2) 

All masses are assumed to be much smaller than the dynamical QCD scale A 

Let us study the extrema of the potential ([T]); our main interest is to follow the dependence 
on the CP- violating phase 9. As we show in the Appendix, one can always search for extrema in 
the diagonal form, U = e'^^^ diag{e'''^^ , e*"^^, . . . e''^^}, with 

N 

^ 0, = + 2Txk , (3) 

i 

where 0j G [— vr, vr) and k is an integer. With this diagonal Ansatz the pion potential becomes 

N 

V{(t)i) = -A^ J2 ^iCos(pi . (4) 



The extremality conditions for the potential (|4]) are simply 

rriiSirKpi = rrijSirKpj, (5) 

for any pair plus the constraint For each value of 9 this set of equations has various 
solutions which correspond to different type of extrema — minima, maxima and saddle points. 

To determine conditions for an extremum to be a local minimum let us consider the Hessian 
matrix Ada/s using 0i, . . . , 4>n-i as independent variables, so that indices a, (3 run from 1 to (A^— 1), 

M-aP = 5ali'maCOS(j)a + rriN COS(f)N, (6) 

where no summation over a is assumed. Then it follows from Eqs. ([2]), ([s]) that a necessary 
condition for a local minimum is cos (pa > 0, a = 1, . . . , A^ — 1. Moreover, if also cos^at > 



^Strictly speaking, in the large Nc regime, one needs <C A/A^c for the description ([l]) to hold, otherwise one 
should also include the rj' meson in the low energy theory. 
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then the extremum is necessary a minimum, while for negative cos (j)N it can be both stable or 
unstable. For equal quark masses rrii = m, this immediately implies that the minima correspond 
to the solutions of (|3]), (|5| where all 0, are equal to each other, 



-e + 27rk 
N 



(7) 



and the parameter k should be in the range 



e 



(8) 
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for an extremum to be a minimum. The value of the axion potential in these minima is 



m cos 



e - 27Tk 
N 



) 



(9) 



At ^ = there are (1 + 2[A^/4]) minima, with k = corresponding to the global one. Let us follow 
now what happens with the energy if one starts at the global minimum at 6' = and adiabatically 
increases the value of 6. As we will see in the next section this is exactly what happens in the 
vicinity of a rotating astrophysical black hole if 9 is promoted to a dynamical axion field with an 
appropriate mass. It follows from ^ that the energy Vq of the k = minimum increases with 9. 
At 9 = 7[ this vacuum ceases to be a global minimum. Finally, there is a critical value 9c, such 
that for 9 > 9c the Vq branch of extrema does not correspond even to a local minimum of energy 
any longer. In the case of equal quark masses 




and Vq is a local maximum for 9 > 9c- In fact, for equal masses it is straightforward to classify 
all the extrema of the pion potential. In Fig. [l^ we reproduced the plot from [22j, where all of 
the extrema are shown for A^ = 3. We see that indeed the branch of extrema corresponding to 
the global minimum at ^ = (for A^ = 3 this is actually the only minimum at ^ = 0) becomes 
metastable at ^ = vr, and turns into a local maximum after hitting a saddle point at ^ = 37r/2. 

A number and properties of extrema change continuously under small variation of the pa- 
rameters of the potential at a generic value of 9. Hence also at non-degenerate masses adiabatic 
variations of 9 may first turn a global minimum into a metastable one, and as the deviation of 
9 from zero increases, may result in a disappearance of a minimum, as soon as the mass split- 
tings are not too large. For instance, for three quark flavors this behavior persists as soon as 
(m2 — m'i)mi < 17121713. Unfortunately, this condition is not satisfied in the real world QCD. 

To illustrate how the structure of extrema changes under variations of masses we presented in 
Fig. different extrema as a function of 9 for mass ratios 1:1.2:1.4. There is still a range of 9 
where a branch corresponding to the global vacuum at ^ = is a local minimum. 

There is one interesting difference of this plot from the one in Fig. [l^. 9 is always a periodic 
parameter with a period equal to 27r. The way this periodicity is realized at equal masses is 
that as 9 varies adiabatically from to 2tt the set of extrema gets reshuffled. This effect is often 
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referred to as axion monodromy. Even at equal masses this name is a bit misleading — as we see 
on the way from to 2tt extrema change their nature and turn from minima into maxima and 
vice versa. Fig. [l^ demonstrates that for non-degenerate masses something even more dramatic 
happens — there are branches of extrema that are not defined over all range of 6 from to 2tt. Of 
course, the physics is still periodic under 27i shifts of 9 — the sets of extrema at 9 and 9 + 2tt are 
identical. 

Note that the two minima have equal energies at 9 = it both for different and degenerate 
masses. This is not a coincidence, but rather a consequence of the 27r-periodicity in 9 and of 
9 — )■ —9 reflection symmetry. 

Let us estimate now the life-time of metastable minima in the pion potential. As usual the 
probability of tunneling is determined by the classical bounce action, F oc e"*^ [231 . It is clear 
that the tunneling probability becomes large for 9 close to 9c, where the height of the barrier goes 
to zero. Indeed the relevant part of the pion Lagrangian at ~ can be written in the form 

£ ~ ^ (50)2 _ {M\9)(f)^ - + . . . ) , (10) 

where M^{9c) = 0, is a linear combination of pions (pi that becomes massless at 9 = 9c, and A 



is a constant. Dots in (10) stand for corrections which are higher order in (j) and/or in {9 — 9c)- 
These terms can be neglected at ^ ~ 6'c when the bounce solution is localized near = C|^ In this 
limit the tunneling rate becomes [25] 

r~exp --^^ . (11) 



^The cubic coupling A can turn zero at some special values of quark masses. In this case one has to include 
the quartic self-interaction to calculate the bounce action. This will change some of our formulas but not the 
conclusions. 
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i.e., the tunneling is indeed unsuppressed for 9 ~ 9c. It will be important for what follows that this 
tunneling regime corresponds to a thick wall limit, when a significant energy ~ mA^ is released in 
pions as a result of a tunneling. For this regime to be physically relevant we need to check that 
the tunneling is slow in the opposite (thin wall) regime. The tunneling probability in the thin 
wall limit is given by [21] 

27 , a' \ 



where AE is the energy difference between the vacua and a is the domain wall tension. To estimate 
a we consider the equal mass case when the transition proceeds between vacua with U = 1 and 
U = e*^'^/^l. The domain wall tension a is then equal to the Euclidean action functional per unit 
surface area 

dx f^Ttd^U^d^U + v\ (13) 

calculated on the optimal solution, which depends on a single variable x. It was calculated for 
three quark flavors in [22]. We will estimate the tension for a general by assuming the same 
Ansatz U{x) = rfmc/{e*^(^), e^^(^) . . . e-^^^-^^^^^)}. The resuh is 



cr 



where In is a numerical coefficient of order 1 (for = 3 it is equal to 6(1 — n / (3^/3)) . This results 
in the following estimate for the tunneling rate. 



16 mA^ sin^ {§) {9 



We see that the thin wall rate is naturally suppressed both by a large numerical factor in the 
exponent, and also because the ratio -FV/m is large in a chiral limit. We conclude that the 
transition should be close to the thick wall regime in order to proceed efficiently, giving rise to a 
significant latent heat ~ mA^, which is released in the form of pions. 

To finish this section note that a non-trivial vacuum structure as a function of 9 parameter of 
the type discussed here appears to be a generic property of strongly coupled gauge theories, not 
limited to a specific calculable example that we chose. In particular, it is likely to be a feature of 
the pure gluodynamics as well 



3 Hidden Valley Avalanches 

Let us describe now how axionic superradiant instability of astrophysical black holes may trigger 
rapid instabilities (avalanches) in hidden valleys. Eventually, our main focus will be on possible 
electromagnetic echoes of these avalanches, but in this section we describe the relevant dynamics 
entirely within a hidden sector. Let us start with a quick review of the basics of superradiant 
instability following Ref. [10]. 
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Light axion fields render rotating astrophysical black holes unstable. The instability rate 
depends on the ratio a of the black hole size and the axion Compton wavelength, a = fia^g- 
For a maximally spinning black hole the rate is maximum, F ~ 10~^r~^, at a ~ 1/2, and rapidly 
decreases when a deviates from this optimal value. In practice the size of astrophysical black 
holes is tiny compared to the age of the Universe, so that a can be significantly different from its 
optimal value and still the instability has enough time to develop. Every axion may affect black 
holes within two-three orders of magnitude in mass (depending on /i^), as illustrated in Fig. 2 of 
Ref. dU]. 

It is useful to think about black hole/axion system as a huge gravitational atom. One crucial 
difference with a hydrogen atom is that axions are bosons so that energy levels may have large 
occupation numbers (of order 10™ in the present context). Another important difference is that 
an atomic nucleus is replaced by a black hole, as a result the boundary condition at the origin 
gets modified, and the energy levels acquire imaginary parts. For a Schwarzschild black hole all 
imaginary parts are negative, implying that a black hole absorbs particles from all levels. For a 
rotating black hole some of the levels acquire positive imaginary parts, signaling the presence of a 
superradiant instability. The occupation numbers of these superradiant levels grow exponentially 
by extracting rotational energy from the ergosphere. This is a variation of the famous Penrose 
process. 

To describe the development of the instability let us start by considering a black hole which 
rotates rapidly at the initial moment of time. This may be a realistic situation if, e.g., a black 
hole is just created as a result of a supernova explosion with a large natal spin. 

From that moment on the black hole starts rapidly loosing its spin by populating the fastest 
superradiant level. The process continues until non-linearities in the axion potential become 
important. For an axion potential of the form shown in Eq. ^ (with k = 0) this happens at 
9/N 1. In the absence of a hidden gauge sector the estimates of Ref. suggest that as an 
outcome of non-linear dynamics the axion cloud collapses — a similar process was observed for 
attractive Bose-Einstein condensates in a lab, and is called "Bosenova" [29]. After the collapse, 
presumably, an order one fraction of axions gets absorbed back into the black hole, and an order one 
fraction escapes. After the environment cleans up the whole process repeats. This happens until 
the black hole spin drops down to the critical value, where the width of the fastest superradiant 
level changes its sign. 

At this point one may expect that the next-to-fastest superradiant level becomes dynamically 
important and the whole story repeats with a slower time-scale set by its population rate. However, 
the presence of the axion cloud populating the first level induces mixing between other levels that 
shuts down the instability for these levels. As a result the process repeats for the second level only 
after a cosmologically long time required for the axion cloud to dissipate, so that the level mixing 
becomes sufficiently small. In the absence of external perturbations the dissipation happens due 
to axion annihilations into gravitons. Alternatively, the cloud can be destabilized as a consequence 
of a merger event, or if the black hole mass and spin changes by order one as a result of a long 
period of accretion. 

In all these cases a time interval between episodes of superradiance instability is cosmologically 
long, so that every black hole may experience several such episodes through the history of the 
Universe. Interestingly, these episodes may be happening for isolated black holes that do not 
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correspond to any visible astronomical sources before the superradiance starts. The focus of 
Ref. [TOj, was mainly on possible gravitational wave signals during the episodes of superradiance 
for a generic axiverse axions and also on the possible non-gravitational signals for the QCD axion. 
Instead, here we consider what happens if an axion causing the superradiance is coupled to the 
QCD-like hidden sector with the phase structure as described in section |2j For concreteness, we 
consider the case of equal quark masses. It is straightforward to extend the whole discussion to 
the case of different but close enough masses, such that the non-trivial phase structure persists. 

In the field theoretic language the growth of the axion cloud implies that the value of the axion 
field gets excited from the minimum of its potential in the vicinity of a black hole. It experiences 
oscillations with a period set by the axion mass, and with an exponentially growing amplitude. 
The spatial profile of the field follows the shape of the hydrogen wave function with a typical 
variation scale set by the Bohr radius, Rc ~ P/{a^a)i where I is the angular momentum of the 
corresponding superradiant level. 

Note that the scales for both time and space variations of the axion field are extremely slow 
compared to the dynamical scale A of the corresponding strongly coupled gauge sector. Conse- 
quently, from the hidden valley viewpoint the interior of the axion cloud can be thought of as 
the QCD laboratory with an adiabatically varying 6'-parameter; the amplitude of this variation 
exponentially (but still extremely slowly) grows in time. 

Eventually, as 9/N becomes larger than n/N, an observer in the cloud finds herself in a 
metastable reheated state. Initially the tunneling rate into the lower vacuum is much longer than 
the axion oscillation time (and, in fact, longer than the age of the Universe), for instance, for 9/N 



close to tt/A^ it is given by the thin wall expression (15). Oscillations continue in this regime until 



their amplitude start approaching 6/N = 7r/2, where mesons become tachyonic. Close to this 



saddle point the nucleation rate is determined by the thick wall result (11) and becomes fast as 6 
approaches the saddle. When the nucleation rate becomes fast enough for a bubble of energetically 
favored vacuum to be created during one period of axion oscillations, the whole overheated part 
of the cloud with 9/N > ir/N experiences the first order phase transition. 

Inside the transition region the latent heat of the overheated QCD vacuum is released in the 
form of mesons. Note that the hidden sector with a dynamical axion has a single vacuum at ^ = 0, 
so the net result of this process is that the axion cloud "boils up" , and creates a cloud of strongly 
interacting mesons inside. Initially at the boundary of this cloud one has a "D-wall" interpolating 
between two QCD vacua with 6 = tt (see, e.g., |28] for a review). In the thick wall tunneling 
regime the energy carried by this wall is subleading with respect to the meson energy, and its 
presence is irrelevant for our purposes. The total energy released in mesons can be estimated as 

f2 

Emes ~ mA'rl ~ j^Mbh , (16) 

where mA^ is the latent heat of the phase transitions, approximately equal to the energy density 
in axions prior to the transition, and Mbh is the black hole mass. Here we approximated the size 
of the cloud by Vg, which corresponds to setting a, I ~ 1. 
Mesons rapidly thermalize at a temperature 

T ~ {mA'y/^ , (17) 
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which is always higher than the meson mass 171^ ~ (mA)^/^, so one ends up generating a gas of 
relativistic (or mildly non-relativistic) strongly interacting particles. 

Clearly, there are significant uncertainties in the estimates (16), (17) but they clearly show 
that for high fa ~ Mqut a significant (10~^ -r- 10^^) fraction of a black hole mass can be converted 
into mesons by this mechanism. Unlike axions, whose couplings are suppressed by the scale fa 
and extremely weak, hidden sector mesons may have appreciable interactions with the Standard 
Model fields, so this energy may be efficiently converted into photons. Before discussing more 
details of this conversion let us comment on one possible caveat behind our reasoning. 

The subtlety is that to trigger the instability the axion field needs to approach large enough 
values 6/N ~ 7r/2. However, as discussed above, the main limiting factor for the growth of the 
axion cloud is the Bosenova collapse that happens at ^/A^ ~ 1 Consequently, one may worry 
that the cloud collapse will always happen before the phase transition in the gauge sector, so that 
black holes can never trigger hidden valley avalanches. 

We believe this pessimistic conclusion is unlikely to be correct. First, even if the Bosenova 
collapse happens before the avalanche, it is very likely that the axion field significantly increases 
in the course of the collapse, so that the avalanche will still happen during the Bosenova event 
itself. To definitely confirm this expectation and to calculate the resulting energy release in pions 
a numerical simulation of the whole process is required, which is beyond the scope of our paper. 
So, as a proof of principle, let us describe an example of a situation where the avalanche definitely 
happens before the Bosenova collapse. 

To achieve this let us add two more parameters to our setup. First, we assume that the axion 
is coupled to the hidden valley topological density with an integer coefficient Na> 1-, 



N f 

Se = J d'x (Pae^-^PTiG^Mx, 



Second, let us assume that the axion potential is dominated by a string contribution rather than 
by the hidden valley. Then the axion potential takes the following form, 

V = A% cos{e - 60) - A^m cos {NaO/N) , (18) 

with A^j > A^m. Then the avalanche transition happens at NaO/N ~ 7r/2, while axion non- 
linearities become important at {9 — 9q) ~ 1. Clearly, there are regions in this extended set of 
parameter, such that the avalanche happens in the regime when the axion field is still linear. Note 
that in this case the energy released in mesons during the avalanche is suppressed by the ratio 
(A^m/A^J as compared to the estimate (16), so this ratio should not be too small. 

Admittedly, this situation may appear more tuned compared to the one we were considering 
to start with. However, as we already stressed, we presented this specific example just as a proof 
of principle, and we expect that avalanches happen under more relaxed conditions as well. 



4 Electromagnetic Echoes of Hidden Valley Avalanches 

Let us turn now to discussing how the hidden valley mesons may be coupled to the Standard 
Model fields and what are the resulting electromagnetic signals in a different range of parameters. 

•^In fact, according to estimates of [10] it is likely to happen even earlier at small a/l — at 6/N ^ 2a/l. 
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There are several possibilities ( "portals" ) of how a hidden valley can communicate to the Standard 
Model fields (see, e.g., [6]). For definiteness, we focus here on the vector portal, but the principal 
qualitative conclusions are model independent. In this case a subgroup of the fiavor group of the 
hidden QCD is gauged, so that some of the hidden sector mesons carry a charge under a new U{1) 
boson A', just like it happens for ordinary mesons. 

In principle, one may consider two different scenario of how A' boson interacts with the Stan- 
dard Model fields. The first possibility is that some of the Standard Model fermions are charged 
under the A', as considered, for instance, in |T]. In this case the collider bounds constrain A' to 
be relatively heavy, mA'/g' ^ 10 TeV. 

The more interesting possibility for our purposes is that the only interaction between the 
Standard Model and A' is due to kinetic mixing with the hypercharge, 

where F's are the hypercharge and A' field strengths and e is the mixing parameter. In this case 
A' can be light or even massless if the mixing parameter e is sufficiently small. One theoretically 
motivated range for mixing is e ~ lO^*^^^^-*. Mixing parameters in this range would be generated at 
one loop by integrating out heavy fields charged under both A' and the hypercharge [30] . However, 
mixing can naturally be much weaker as well. This would happen, for instance, if A' is a part 
of a compact group, which is spontaneously broken at a scale much lower than the masses of all 
fields carrying both A' and Standard Model charges. We will treat e as a free phenomenological 
parameter. The summary of observational constraints on e for different range of A' and hidden 
pion masses can be found in, e.g., |8j. 

Let us see now how the constraints on the strength of interactions between a hidden sector 
and the Standard Model depend on the mass scale of a hidden gauge group. More precisely, 
the relevant parameter is A^m, which determines the overall scale of the hidden valley potential 



and, as a consequence, the initial temperature of the meson cloud (17), the frequency scale of the 
resulting electromagnetic radiation, and also the value of the axion mass, yU^ ~ A^m//^. Only 
relatively light sectors with A^m < (IGeV)^ are relevant for our purposes, because axion masses 
for hidden valleys with larger energy scales are too heavy to trigger the superradiant instability 
for astrophysical black holes. 

Regardless of a specific communication portal there are two different physical regimes, depend- 
ing on the masses of hidden pions. Namely, pions lighter than a few MeV may contribute as light 
species to the total energy budget during the Big Bang Nucleosynthesis (BBN). If the confinement 
scale A itself is that low, the same problem arises due to hidden gluons and quarks. To preserve 
the success of the BBN a hidden valley needs to be very weakly coupled to the Standard Model 
fields in this case, so that hidden mesons never reach thermal equilibrium after reheating in the 
early Universe. Typically, there is also a strict constraint from above on the reheating temperature 
after infiation, because at high temperatures thermalization through higher dimensional operators 
(for instance, related to heavy A' exchanges) becomes more efficient. 

For a mass scale of a hidden sector higher than ~ 10 MeV the constraints described above get 
significantly relaxed and interactions between the Standard Model and hidden sector fields may 
be much stronger. In particular, it is possible for the hidden sector to be thermalized in the early 
Universe. In this case late decays of the hidden sector fields may cause distortions to the CMB 
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black body spectrum or affect the BBN predictions if the hfe-time is too long. This issue does not 
arise if the life-time is shorter than ~ 1 s. This leads us to consider the scenario with a kinetically 
mixed light A', otherwise the life-time gets too long for sub-GeV hidden valleys which we consider 
[1]. Then one can take e as high as ~ 10~^ -7- 10^^. A particularly interesting region of parameters 
is when the A' mass is also in the range ~ 10 ^ 10'^ MeV, which can also be probed in the near 
future fixed-target experiments [5]. This mass scale of the hidden sector corresponds to stellar 
mass black holes. 

For lighter hidden sector mesons their mixing with the Standard Model should be very sup- 
pressed. The natural option to consider in this case is that the hidden photon is massless, so that 
hidden mesons acquire a tiny electric charge e in the diagonal basis for photons. Observational 
bounds on the mixing parameter e are very strict for light millicharged particles [31118]. For meson 
masses smaller than ~few MeV the mixing parameter should be smaller than ~ 10~^ from the 
cooling of SN1987A, for masses below ~100 keV — smaller than ~ 10~^^ from the white dwarfs 
cooling, and smaller than ~ lO"^'* for masses below ~ 10 keV from red giants cooling. 

Note, however, that for light hidden sectors axions also gets lighter and as a result the size of 
black holes affected by superradiance increases. Consequently the total amount of emitted elec- 
tromagnetic energy can still be very significant. Another reason making this scenario potentially 
interesting is that the detectors sensitivity generally increases when frequencies are getting lower. 

The size of a cloud (in these estimates we do not distinguish between the radii of a black hole 
and of a cloud) scales with the characteristic temperature T ~ (A^m)^/^ as 

where RgQ ^ lA km is the gravitational radius of the Sun. From this estimate we see that the 
most energetic photons originating from pion decays and collisions may have frequencies ranging 
between ~ 100 MeV for stellar mass black holes and ~ 100 eV for the largest supermassive black 
holes with masses of order lO^^M©. 

Depending on the interaction strength and the size of the cloud, two evolution scenarios are 
possible. The distinction is whether hidden sector particles are thermalized with e~^e~ and photons 
or not. The thermalization condition reads 

naR. ~ eV^Tfi. ~ 7 • 10«.^ (^) > 1 • (21) 

If this inequality is satisfied (which is possible for meson masses larger than ~ 100 keV) then the 
cloud behaves as a thermal gas including an electromagnetic component with an initial temper- 
ature T ~ (A^m)^/^ ~ (/i/a)^/^. This gas starts to expand and cool under its own pressure. So 
a likely outcome of the avalanche is a creation of a fireball consisting of photons and e~^e~ pairs 



and carrying according to (16) up to ~ 10""^"^^ of a black hole mass. 

Interestingly, this is similar to the starting point of one of the most successful model (the 
fireball model) for the 7-ray burst (see, e.g., [521 for reviews). The required fireball energy in 
this model is comparable to the avalanche energies in our setup, which is not surprising given that 
both energies are set by a rest mass of a stellar size object. Also the rate of avalanches is likely 
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to be close to the rate of 7-ray bursts. Indeed, a significant fraction of 7-ray bursts is likely to 
originate from supernovae explosions leading to a birth of a rapidly spinning black hole [M]. In 
turn, every spinning black hole is expected to experience several periods of superradiant activity 
per the age of the Universe in the presence of an appropriate axion. 

There is an important difference between our setup and traditional ones, however. An im- 
portant ingredient in the fireball model for 7-ray bursts is the baryonic "contamination" of the 
electromagnetic plasma. This results in a significant delay of the photon signal, because the 
electromagnetic energy first gets converted into the kinetic energy of baryons and only later gets 
released into photons as a consequence of a shock formation. Given that superradiance may 
happen for a single black hole long after the initial supernovae explosion it is possible to have 
practically no baryon contamination in our case. The corresponding avalanches are likely to give 
rise to 7-ray burst-like events of very short duration and with a spectrum close to thermal. 

Another important property of conventional 7-ray burst models is a strong coUimation of 
electromagnetic energy into relativistic jets. The initial shape of a superradiant cloud is highly 
non-spherical, and it is formed in a close vicinity of a rapidly rotating black hole, so it might 
give rise to a jet as well. Also, as we discussed, superradiance is likely to give rise to a series of 
repeating Bosenova collapses, so that these bursts may repeat themselves. 



Let us discuss now what happens if the thermalization condition (21) is not satisfied (which 
must be the case for mesons lighter than ~ 100 keV). In this case after mesons cool down they 
all decay into hidden photons, so the only electromagnetic signal comes from the hot phase when 
ordinary photons are created in the hidden sector plasma by e-suppressed processes. As before, 
the life-time for the hot phase is set by the cloud size Rc 

To estimate the resulting electromagnetic energy release in this case let us denote by A^ the 
maximum energy scale of a hidden sector that can be reheated by an axion with a mass n, 

The actual temperature of mesons in the cloud can be lower than A^ if the dominant contribution 



to the axion potential comes from some other source, as in (18). Then the electromagnetic energy 
release during the time interval ~ Rc can be estimated as 

where we assumed that photons are produced predominantly as a result of bremsstrahlung, and 
estimated the corresponding thermal cross-section as cr ~ e'^aEAiT''^. 

This estimate assumes that the total energy of produced photons is smaller than the total 
energy of the cloud 

Mc ~ T^Rl (23) 



which is indeed the case if the thermalization condition (21) does not hold. 

To illustrate the potentially observable effects of the hidden sector avalanches in Fig. |2] we plot 
an estimated total energy release as a function of the characteristic temperature of radiation. On 
the same plot we show the corresponding black hole masses and the characteristic time scales. 



"^If the hidden photon is massive with a Hfe-time much longer than Rc, there is also a delayed electromagnetic 
signal from hidden photon decays. 
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Figure 2: A total energy release as a function of a characteristic temperature of a fireball for 
different values of parameters a and /. Two red axes on top represent a corresponding black hole 
mass and time scale for a = 0.5, / = 1. For the blue line a black hole mass is 5 times larger than 
the value on the red axis and a time scale is 5 times longer. For the green line a black hole is 5 
times lighter and a time scale is ~ 5^^^ times longer. 



To produce this plot we assumed that the axion potential is dominated by the hidden sector, 
so that T = A„ with fa = 2 ■ 10^® GeV. The total energy release is determined either by (22) or 



by the total mass of the axion cloud (23) depending on whether thermalization condition holds or 
not. To check the latter we assumed that the mixing parameter e takes the largest possible value 
compatible with the observational bounds at the corresponding meson mass scale. 

The reason we have three different curves on the plot is the following. In our previous estimates 
we were assuming a ~ / ~ 1. The three curves are meant to illustrate the uncertainties in the 
estimates related to variations of these two parameters. The red solid curve on the plot corresponds 
to I = 1 a = 1/2, i.e., to a ^ I ^ 1 regim^ The dependence on a and / comes mainly from 
two effects. First, the size of the cloud scales as Rc ~ l^/a, affecting in turn the energy release 



through expressions (22), (23) and through the i?c-dependence of the thermalization condition. 
To illustrate this effect we presented the blue dashed curve, corresponding to Z = 5, a = 2.5. 

Yet another effect enters when we vary a and /, such that the ratio a/l changes as well. 
Namely, as we already mentioned in section |3} the value of the 6 parameter when the Bosenova 
collapse starts was estimated in Ref. [10] to be proportional to a/l. Consequently, for small a/l 



^Recall that the superradiance condition imphes that a/l < 1/2. 
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the avalanche definitely does not happen before the Bosenova. Assuming that it still happens in 
the course of the Bosenova we can estimate the energy release in the following way. Let us take 
as a simplified model of an initial stage of Bosenova that the axion cloud shrinks while preserving 
the total number of axions N ~ O^R^. If this is the case, 6 inside the cloud grows and reaches the 
value 9c when the size of the cloud is R'^ ~ {a/0.5)^/^Rc. Using this size instead of Rc to estimate 
the energy release we obtain the dotted green curve on the plot, which corresponds to / = 1, 
a = 0.1. Clearly, uncertainties for this curve are even stronger than for our previous estimates. 

5 Discussion 

In this paper we outlined a proposal how rotating astrophysical black hole may open a new window 
into hidden QCD-like gauge sectors coupled to light axions. Thanks to axion superradiance rotat- 
ing black holes may give rise to violent bursts of electromagnetic radiation with initial frequencies 
set by the confinement scale in the hidden sector. Clearly, a number of details need to be filled 
in to extract concrete predictions for the characteristics of these sources, allowing to distinguish 
them from astrophysical sources of a more conventional origin and to set limits on the existence 
of these phenomena at different frequencies. 

First, our proposal relies on a qualitative scenario for the evolution of superradiance instability 
developed in Ref. [10]. Given that a black hole surrounded by an axionic cloud is a very rich 
non-linear dynamical system it is highly desirable to confirm this scenario by detailed numerical 
simulations. 

However, many features of our proposal can be studied without going into details of super- 
radiance. Namely, from a purely phenomenological point of view what happens is that a black 
hole creates a fireball of photons and, if the frequencies are high enough, of electron/positron 
pairs. The initial temperature of the fireball is bounded from above by ~ (fa/RgY^'^, where Rg 
is the black hole size and fa is the axion decay constant. This is a rather common situation in 
astrophysics, for instance the initial stage of the fireball model for 7-ray bursts is very similar. 

There is a number of differences in our case which affect the further evolution of the fireball and 
should help to distinguish this kind of sources. For instance, according to the results of Ref. [TU] , 
an episode of superradiant instability may happen for a single black hole in a medium with a very 
small baryon contamination. In this case we expect the resulting burst to have spectrum close to 
thermal and have a shorter time scale than a typical 7-ray burst (here we refer to the case, when 
the frequency scale of the event corresponds to 7-rays). On the other hand, the electromagnetic 
fireball suggested here is immersed in a cloud of a dense liquid composed of strongly interacting 
hidden mesons, which may affect its evolution. 

When the episode of superradiance happens to a black hole which was just created as a result 
of supernova explosion, the whole setup becomes really similar to conventional 7-burst models. 
Even in this case the details may be different, though. For instance, superradiance is a relatively 
slow process. Furthermore, it can be shut off in a dense medium, so that it may be necessary 
for an environment to clean up a bit for the instability to start. As a consequence we expect a 
significant time delay between a black hole collapse and the creation of a fireball. 

Finally, supermassive black holes may give rise through superradiance to very powerful elec- 
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tromagnetic bursts at much lower (down to ~ 100 eV) frequencies than conventional 7-ray bursts. 

To conclude, superradiance opens an intriguing opportunity to test new physics with astro- 
physical black holes. Given that violent electromagnetic bursts are being observed in nature and 
the origin of some of them is unknown we beheve a dedicated astrophysical study of the process 
proposed here is worthwhile. 
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Appendix: Extrema of the pion potenial are flavor diagonal 

Let us first show that for non-degenerate quark masses all the extrema correspond to diagonal U 
matrices. By taking the derivative of ([T]) we arrive at the following extremality conditions, 

Tr [i {e-''/^MU - U^e^'/^M) Ha) = (24) 

where Ha form an orthonormal basis in the space of traceless Hermitian matrices. This condition 
implies that 

e-''/^MU - f/te^'/^M = irl, (25) 



for some real number r. Taking the square of (25) and substituting 



which also follows from ( 25 ) , we get 

UqM^ = M'^Uq. (26) 

Clearly, if all elements of a diagonal matrix M are different, a unitary matrix also has to be 
diagonal to commute with it. 

Finally, if k of the quark masses are degenerate, an extremum may have a non-diagonal k x k 
block. However, this block can be diagonalized by a vector SU{k)v rotation of mesons, which is a 
symmetry of the chiral Lagrangian in this case. Consequently, we can always search for extrema 
in the diagonal form. 



References 

[1] M. J. Strassler and K. M. Zurek, Phys. Lett. B 651, 374 (2007) [arXiv:hep-ph/0604261]. 

[2] A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper and J. March-Russell, Phys. Rev. 
D81, 123530 (2010). [arXiv:0905.4720 [hep-th]]. 



15 



[3] J. Kang, M. A. Luty, JHEP 0911, 065 (2009). [arXiv:0805.4642 [hep-ph]]. 

[4] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer and N. Weiner Phys. Rev. D79, 015014 
(2009). [arXiv:0810.0713 [hep-ph]]. 

[5] J. D. Bjorken, R. Essig, P. Schuster and N. Toro, Phys. Rev. D 80, 075018 (2009) 
[arXiv:0906.0580 [hep-ph]]. 

[6] B. Batell, M. Pospelov and A. Ritz, Phys. Rev. D 80, 095024 (2009) [arXiv:0906.5614 [hep- 
ph]]. 

[7] A. Arvanitaki, N. Craig, S. Dimopoulos, S. Dubovsky and J. March-RusseU, Phys. Rev. D81, 
075018 (2010). [arXiv:0909.5440 [hep-ph]]. 

[8] J. Jaeckel, A. Ringwald, Ann. Rev. NucL Part. Sci. 60, 405-437 (2010). [arXiv: 1002.0329 
[hep-ph]]. 

[9] C. Cheung, Y. Nomura, J. Thaler, JHEP 1003, 073 (2010). [arXiv: 1002. 1967 [hep-ph]]. 

[10] A. Arvanitaki and S. Dubovsky, arXiv: 1004.3558 [hep-th]. 

[11] R. Penrose, Riv. Nuovo Cim. 1, 252 (1969) [Gen. Rel. Grav. 34, 1141 (2002)]. 

[12] D. Christodoulou, Phys. Rev. Lett. 25, 1596 (1970). 

[13] Ya. B. Zeldovich, JETP Lett. 14, 180 (1971) 

[14] C. W. Misner, Phys. Rev. Lett. 28, 994 (1972). 

[15] A. A. Starobinskii, Soviet Physics JETP, 37, 28 (1973) 

[16] T. Damour, N. Deruelle and R. Ruffini, Lett. Nuovo Cim. 15, 257 (1976). 

[17] I. M. Ternov, V. R. KhaUlov, G. A. Chizhov and A. B. Gaina, Sov. Phys. J. 21, 1200 (1978) 
[Izv. Vuz. Fiz. 21N9, 109 (1978)]. 

[18] T. J. M. Zouros and D. M. Eardley, Annals Phys. 118, 139 (1979). 

[19] S. Detweiler, Phys. Rev. D 22, 2323 (1980). 

[20] E. Witten, Annals Phys. 128, 363 (1980). 

[21] M. Creutz, Phys. Rev. D52, 2951-2959 (1995). [hep-th/9505112]. 

[22] A. V. Smilga, Phys. Rev. D59, 114021 (1999). [hep-ph/9805214]. 

[23] I. Y. Kobzarev, L. B. Okun, M. B. Voloshin, Sov. J. Nucl. Phys. 20, 644-646 (1975). 

[24] S. R. Coleman, Phys. Rev. D15, 2929-2936 (1977). 



16 



[25] A. D. Linde, Nucl. Phys. B216, 421 (1983). 

[26] E. Witten, Phys. Rev. Lett. 81, 2862-2865 (1998). [hep-th/9807109]. 
[27] M. A. Shifman, Phys. Rev. D59, 021501 (1999). [hep-th/9809184]. 

[28] G. Gabadadze, M. Shifman, Int. J. Mod. Phys. A17, 3689-3728 (2002). [hep-ph/0206123]. 

[29] E. A. Donley, N. R. Claussen, S. L. Cornish, J. L. Roberts, E. A. Cornell, and C. E. Wieman, 
Nature 412, 295(2001) 

[30] B. Holdom, Phys. Lett. B178, 65 (1986). 

[31] S. Davidson, S. Hannestad, G. Raffelt, JHEP 0005, 003 (2000). [hep-ph/0001179]. 

[32] T. Piran, Rev. Mod. Phys. 76, 1143-1210 (2004). [astro-ph/0405503]. 

[33] P. Meszaros, Rept. Prog. Phys. 69, 2259-2322 (2006). [astro-ph/0605208]. 

[34] A. MacFadyen, S. E. Woosley, Astrophys. J. 524, 262 (1999). [astro-ph/9810274]. 



17 



